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We study decoherence induced by stochastic squeezing control errors considering the particular 
implementation of Hadamard gate on optical and ion trap holonomic quantum computers. We 
analytically obtain both the purity of the final state and the fidelity for Hadamard gate when the 
control noise is modeled by Ornstein-Uhlenbeck stochastic process. We demonstrate the purity and 
the fidelity oscillations depending on the choice of the initial superimposed state. We derive a linear 
formulae connecting the gate fidelity and the purity of the final state. 
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I. INTRODUCTION 

Holonomic quantum computations exploiting non- 
abelian geometrical phases [jj was primarily proposed 
in the Ref. 2] and developed further in the Ref. 0. 
Many implementations of holonomic quantum comput- 
ers (HQC) have been proposed. Particularly, the real- 
ization of HQC within quantum optics was suggested 
(optical HQC) 4]. Laser beams in a non-linear Kerr 
medium were exploited for this purpose. Two different 
sets of control devices can be used in this case. The 
first one considered in this paper consists of one- and 
two-mode displacing and squeezing devices. The second 
one includes SU(2) interferometers. As well trapped ions 
with the excited state connected to a triple degenerate 
subspace (four level A-system) can be used to imple- 
ment HQC 5]. Another approach to HQC exploiting 
squeezing and displacement of the trapped ions vibra- 
tional modes was suggested in the Ref. 6|. This im- 
plementation of HQC is mathematically similar to the 
first embodiment of the optical HQC Q| and thus it is 
also considered in this work. Particularly, expressions for 
the adiabatic connection and holonomies are the same in 
these cases. Another proposed implementation of HQC 
was the HQC with neutral atoms in cavity QED . The 
coding space was spanned by the dark states of the atom 
trapped in a cavity. Dynamics of the atom was gov- 
erned by the generalized A-system Hamiltonian. Math- 
ematically similar semiconductor-based implementation 
of HQC was proposed in the Refs. @, where one-qubit 
gates were also realized in the framework of the general- 
ized A-system. In distinction from the cavity model of 
HQC its physical implementation exploits semiconduc- 
tor excitons driven by sequences of laser pulses . For 
the two-qubit gate implementation the bi-excitonic shift 
was used. The generalized A-system with the different 
Rabi frequencies parametrization was exploited recently 
for HQC implemented by Rf-SQUIDs coupled through a 



microwave cavity 0. One more solid state implementa- 
tion of HQC based on Stark effect was proposed in the 
Ref. [13. 

Let us briefly remind the main results concerning the 
holonomic quantum computation. In HQC non-abelian 
geometric phases (holonomies) are exploited to imple- 
ment unitary transformations over the quantum code. 
The later is some degenerate subspace C N spanned on 
eigenvectors of Hamiltonian Hq, which initiates the para- 
metric isospectral family of Hamiltonians F — {H(\) = 
U(X)HoU'(\)}\ e M- Here U(X) is a unitary operator, A 
is a vector belonging to the space of the control param- 
eters M and N denotes the dimension of the degenerate 
computational subspace 0, Q . Quantum gates are im- 
plemented when the control parameters are adiabatically 
driven along the loops in the control manifold M. The 
unitary operator mapping the initial state vector belong- 
ing to C N into the final one has the form e l ^T 7 (A^), 
where the index fi enumerates control parameters, A^ 
constitute vector A and <fi is the dynamical phase. Holon- 
omy associated with the loop 7 S M is 



r 7 (^) = Pexp 



(1) 



P denotes the path ordering operator, is the 



Here 

matrix valued adiabatic connection given by the expres- 
sion u 



{<Pr, 



(2) 
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where \^k) with k — 1,N are the eigenvectors of the 
Hamiltonian H forming the basis in C N . Dynamical 
phase <f> will be omitted bellow due to the suitable choice 
of the zero energy level. We shall consider the single 
subspace C N (no energy level crossings are assumed). 

It is evident that the quantum gate (holonomy) per- 
formed depends on the path passed in the control param- 
eters space. As well it is obvious that in real experiments 
it is impossible to pass the desired loop in the control 
manifold without any deviations. Errors in the assign- 
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ment of the classical control parameters A are unavoid- 
able. The question about robustness of holonomic quan- 
tum computations with respect to the control errors has 
attracted a lot of attention recently. Namely, the effect of 
the errors originated from the imperfect control of clas- 
sical parameters was studied for CP™ model of HQC in 
the Ref. where the control-not and Hadamard gates 
were particularly considered. Berry phase for the spin 
1/2 particle in a classical fluctuating magnetic field was 
considered in the Ref. 1131. Approach based on the non- 
abelian Stokes theorem |l3j] was proposed in the Ref. [l4| . 
Namely, the general expression for the fidelity valid for 
arbitrary implementation of HQC in the case of the single 
control error having arbitrary size and duration was de- 
rived. Simple approximate formulae was obtained in the 
small error limit. Adiabatic dynamics of quantum system 
coupled to a noisy classical control field was studied in the 
Ref. . It was demonstrated that stochastic phase shift 
arising in the off-diagonal elements of the system's den- 
sity matrix can cause decoherence. The efficiency of Shor 
algorithm [HI run on a geometric quantum computer was 
investigated in the case when the decoherence induced 
by the stochastic control errors was taken into account. 
The study of the robustness of the non-abelian holonomic 
quantum gates with respect to the stochastic fluctuations 
of the control parameters was presented in the Ref. [l7| . 
Three stability regimes were discriminated in this work 
for the HQC model with qubits given by polarized exci- 
tonic states controlled by laser pulses. Noise cancellation 
effect for simple quantum systems was considered in the 
Ref. 0|- Robustness of the parametric family of quan- 
tum gates subjected to stochastic fluctuations of the con- 
trol parameters was studied in the Ref. 0] . Usage of the 
cyclic states allowed to consider quantum gates which 
could be continuously changed from dynamic gates to 
purely geometric ones. It was shown that the maximum 
of the gate fidelity corresponds to quantum gates with a 
vanishing dynamical phase. Robust Hadamard gate im- 
plementation for optical Q and ion trap holonomic 
quantum computers was proposed in the Ref. plj . The 
cancellation of the small squeezing control errors up to 
the fourth order on their magnitude was demonstrated. 
Hadamard gate is one of the key elements of the main 
quantum algorithms, for instance see [iH l22|. Thus the 
search for its robust implementations is of importance. 

During the last few years much attention has been 
payed to the study of both abelian and non-abelian ge- 
ometric phases in the presence of decoherence which is 
the most important limiting factor for quantum compu- 
tations. Let us briefly overview some of these works. 
The abelian geometric phase of the two-level quantum 
system interacting with a one and two mode quantum 
field subjected to the decoherence was considered in the 
Ref. [23|. It was demonstrated that when the geomet- 
ric phase is generated by an adiabatic evolution the first 
correction due to the decoherence of the driving quan- 
tized field for the no-jump trajectory has the second or- 
der in the decaying rate of the field but it is not the case 



for the non-adiabatic evolution. Non-abelian holonomies 
in the presence of decoherence were investigated in the 
Ref. [24| using the quantum jump approach. The effects 
of environment on a universal set of holonomic quantum 
gates were analyzed. Refocusing schemes for holonomic 
quantum computation in the presence of dissipation were 
discussed in the Ref. (2{|. It has been shown that non- 
abelian geometric gates realized by means of refocused 
double-loop scheme possessed a certain resilience against 
decoherence. Quantum Langevin approach has been used 
to study the evolution of two-level system with a slowly 
varying Hamiltonian and interacting with a quantum en- 
vironment modeled as a bath of harmonic oscillators (2(| . 
It allowed to obtain the dissipation time and the correc- 
tion to Berry phase in the case of adiabatic cyclic evolu- 
tion. The realization of universal set of holonomic quan- 
tum gates acting on decoherence-free subspaces has been 
proposed in the Ref. j2?|]. It has been shown how it can 
be implemented in the contexts of trapped ions and quan- 
tum dots. The performance of holonomic quantum gates 
in semi-conductor quantum dots under the effect of dis- 
sipative environment has been studied in the Ref. [2£| . It 
was demonstrated that the influence of the environment 
modeled by the superhomic thermal bath of harmonic 
oscillators could be practically suppressed. The study of 
the non-adiabatic dynamics and effects of quantum noise 
for the ion trap setup proposed in the Ref. [5j has been 
also done p9| . The optimal finite operation time was de- 
termined. In the references mentioned above the fidelity 
was used as the main measure of gate resilience. 

In this paper we consider optical and ion trap imple- 
mentations of HQC proposed in the Refs. £| and 
respectively. Regarding the particular implementation 
of Hadamard gate we study the decoherence induced 
by stochastic squeezing control errors. Following the 
Ref. ^2 we model the random fluctuations by Ornstein- 
Uhlenbeck stochastic process. We analytically obtain the 
final state purity and the gate fidelity as the measures of 
the gate robustness with respect to the decoherence in- 
duced by stochastic control errors. In the small squeezing 
control errors limit we derive a simple formulae connect- 
ing the purity of the final state and the gate fidelity. 

II. HADAMARD GATE IMPLEMENTATION 

Optical and ion trap setups of HQC are mathemati- 
cally equivalent since the corresponding holonomies are 
the same (compare the Refs. and Therefore we 
can consider both HQC models simultaneously. 

The laser beams in the nonlinear Kerr medium are 
explored in order to perform holonomic quantum com- 
putation in the framework of the optical setup. The cor- 
responding interaction Hamiltonian describing a single 
beam in the medium is 

= hXata (a f a - l) , (3) 

where a and a* are the annihilation and creation opera- 
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tors of the photons respectively, X is a constant propor- 
tional to the third order nonlinear susceptibility of the 
medium. The degenerate computational subspace of the 
single qubit is spanned on the photon Fock states |0) and 
|1). More details one can find in the Ref. p|. 

In order to implement the same holonomic quantum 
computational scheme in the framework of the ion trap 
setup one has to deal with the two-level trapped ion 
placed in the common node of two standing electromag- 
netic waves with frequencies loq~ oj z and ujq+uj z as well as 
being affected by the traveling wave with the frequency 
loq. Here loq denotes the frequency corresponding to the 
transition between the two ion levels being exploited, w z 
is the frequency of the ion's harmonic oscillations along 
the z axis of the linear Paul trap. The basis qubit states 
are \g) <g> |0) and (\g) ® |1) - |e) <g> |0)) /V2. Here \g) and 
|e) are the ground and excited internal states of the ion, 
|0) and |1) are the two lowest vibrational Fock states of 
the ion in the trap. More details can be found in the 
Ref. 0. 

One-qubit gates are given as a sequence of single mode 
squeezing and displacing operations 0, : 



U(t],u) =D(r,)S(v), 



where 



,t2 



S(y) — exp [va'~ — v a~ 
D(rj) = exp (rja) — rfa) 



(4) 



(5) 



denote single mode squeezing and displacing operators 
respectively, v — r\e tBl - and r\ = x + iy are correspond- 
ing complex control parameters, a and are annihilation 
and creation operators. The asterix denotes complex con- 
jugation. The expressions for the adiabatic connection 
and the curvature tensor can be found in the Refs. [H|(|. 
Following our previous Letter pH we consider Hadamard 
gate 



H 



1 

V2 



1 1 

1 -1 



(6) 



implemented when two rectangular loops belonging to 
the planes (x, ri)|e 1= o and (y, ri) |e 1= o are passed. 
Namely, 

— iH = r(Cjj) |sj/=7r/2 r(C/) |£ /=7r /4 , (7) 

where the holonomies are 

r(d) =exp(-ia y Zi), Sj= J dxdr x 2e-' lr \ 

S(C) 

T(C n ) = exp (-ia^n), Su= J dydr x 2e 2r \ (8) 

S{C n ) 

and S(Ciji) are the regions in the planes {x,r\) le^o 
and (y, r\) |g 1= o enclosed by the rectangular loops Cj and 
Cu respectively. The sides of the rectangles Cj and Cu 



are parallel to the coordinate axes. For the loop C/ these 
sides are given by the lines r\ = 0, x — b Xl r\ — d x and 
x = a x , where the length of the rectangle's sides parallel 
to the x axis is l x = b x — a x . In the Ref. |2lj it was shown 
that 



1 / 7T 



I x > -■ 

4 



(9) 



In the same way the rectangle Cji is composed of the 
lines r\ = 0, y = b y , r\ = d y and y = a y , where |2lj : 



1, /, 7T 

(L = - In 1H 

y 2 \ 2L 



ly by dy . 



(10) 



Proposed Hadamard gate implementation is not a unique 
one. The same gate can be realized by passing another 
loops in the control manifold. Our choice is motivated 
by the simplicity of the loops. 



III. DECOHERENCE INDUCED BY 
STOCHASTIC SQUEEZING CONTROL ERRORS 

We restrict ourselves by the consideration of the 
squeezing control errors only. Moreover, we can neglect 
the fluctuations of the squeezing control parameter when 
7*1 = 0. Thus to take into account random squeezing con- 
trol errors we have to replace d x by d x + Sr x (x) and d y 
by d y + Sr y (y), where 5r x (x) and Sr y (y) are two indepen- 
dent Ornstein-Uhlenbeck stochastic processes. Making 
this substitution into the Eqs. JSJ instead of the formu- 
lae (JZJ we obtain the following expression for the per- 
turbed Hadamard gate, see also |21j : 



— iH = 7= (cos a — sin a) (sin/3 + ia x cos/3) — 

v2 

1 

-= (cos a + sin a) (a z cos (3 — a y sin 0) , (11) 

where 



a = e 2dx dx (l — e 



13 = e 2d y / dy (e 2Ar « - l) 



(12) 



Let the qubit initially to be in the pure superimposed 
state |-0o) = c o |0) + C\ |1), where amplitudes Cq and Ci 
obey the normalization constrain |co| 2 + |ci| 2 = 1. For 
the fixed noise realization the final qubit state will be 
pure. However, it will differ from the desired one. In the 
real experiment we do not follow the random fluctua- 
tions of the control parameters (nevertheless in principle 
we can do it). In this situation quantum mechanics pre- 
scribes us to describe the final state of the system by the 
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density matrix and represent the state as a mixture of 
all possible final states weighted with the probabilities 
of the corresponding noise realizations. Following this 
strategy we find the density matrix of the final state for 
a given noise implementation and than average over the 



squeezing control parameter fluctuations when the later 
are modeled by the two independent Ornstein-Uhlenbeck 
stochastic processes. 

Thus, for the density operator p = H \ipo) (ipo \ we 
obtain the following matrix elements: 



(°M°) = 2 + 2 U C °' 2 ~~ ' Cl cos2 P cos2 7 ~ 2 ( c o c i +CoCi)cos2/3sin27- - (cqc\ - cjci) sin 2/3 ', 

l/2 2\ i f 2 2\ 1 

<°l Pi 1 ) = 2 V Cl l — I c o | J sin27 + - [\c \ - \a\ J sin2/3cos27 - - (c c\ + c Ci) (cos27 + i sin 2/3 sin 27) 

~\ ( c oc* ~ c* Cl )cos2/3, 
(l|p|0) = (0|p|ir, (l|p|l) = l-(0|p|0). (13) 



Here the introduced parameter 7 is defined as 7 = a— 7r/4 
and the asterix denotes the complex conjugate quantities. 

We assume that the noise Sr x has variance a x and a 
lorentzian spectrum with the bandwidth T x . The fluc- 
tuations Sr y have the variance a y and bandwidth T y . 
Using the Eqs. ifT^ - lTH?)) and the properties of Ornstein- 
Uhlenbeck stochastic process (see Ref. pfflj]!: 

Sr Xt y(h)Sr x ^{t 2 ) = c^e-^l* 1 -* 2 !, 

Srl :y = a- Xt y, Sr Xt y = 0, (14) 



we average the density matrix p over the stochastic fluc- 
tuations of the squeezing control parameters Sr x and 6r y . 
The line over the random quantities means the averag- 
ing operation. We assume that Sr x ^ y <C 1 and restricted 
ourselves by the first non-vanishing terms depending on 
5r x or 5r y . As the result of straightforward but a bit 
lengthy calculations we analytically obtain the elements 
of the averaged density matrix ~p: 



I P |0> = ^ (1 + coc*, + cgci) - 2 ( I c 1 2 - |ci| 2 ) aj x 



2i (c c* - c* Cl ) a y l y e 2d y 



(cqC* + CqCi) 



8a x 



-T 



I P |1> = 2 ( c o + c i) ( c o - C D + 2 (co c i + c&ci) (^^e" 2 ^ + ia y l y e 2d ») 

\ci\ 2 ) T x + ^ (cqcI - clci)T Vl 

I 1 y 



— (M 

J- a; x 

(l|p|0) = (0|p|l)% <1| p |i> = 1 - <0| p|0> 

I 



(15) 



Here we introduced the following denotations: 

1 _ e -r.i. • 



-4d x 



r, 

1 - e 



— r 1 



r„ 



(16) 



In order to quantify decoherence strength we exploit 
the purity of the final state. It is defined as the trace of 
the squared density matrix. Purity equals to 1 for pure 
states and less than 1 overwise. Using the Eqs. (|13fl it is 
easy to check that for a fixed noise realization the purity 



Iq = trp 2 equals the unity as it should be. Using the 
Eqs. (|15|) we obtain the purity of the final state in the 
case of the stochastic squeezing control errors. The result 
of the lengthy but straightforward calculations is 



: trp 2 
I6a x 



1 



64er t 



Fy\c \ \ci\ 



Fx (Cn 



(17) 



Thus we see that the final state purity I < 1 and the 
stochastic squeezing control errors induce decoherence 
and lead the final state to be a mixture of pure states. 
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We can simplify the expression 1|17|) if we exploit the fol- 
lowing parametrization for the amplitudes Cq and c±: 



cq = e 1 ^ cos if, c\ — e* x sin ip, 



(18) 



and assume that £ — x — nn , where n is an integer. In 
this rather general case formulae l|17|) reduces to 



/ = 1 



16<7 T 16(7,, o 

— -T x - -^T y sin 2 2y. (19) 

^ X L y 



We see that the final state purity oscillates depending on 
the choice of the initial qubit state. The purity decreases 
when the initial qubit state is a superposition of the basis 



states. For example, it has its minimum when \ipo) = 
(|0) ± |1)) j\[2. In contrary, the purity has its maximum 
if the initial qubit state |^>o) is proportional to one of the 
basis states |0) or |1). In this case the errors made in the 
{y-i r i) |fli=o plane are eliminated. 



IV. HADAMARD GATE FIDELITY 

Now we obtain the fidelity of the non-ideal Hadamard 
gate. In the case when there are no control errors (Sr x = 
Sr y — 0) the density matrix p 1 - ^ of the final (pure) state 



has the following form: 



J 



c ci 

|2 



I CO 



|co| - \c%\ +cqc\-cIci 



c ci 



c ci - coCi 
- c cl 



(20) 



r 



The non-ideal Hadamard gate fidelity F = tr(p(°'p) un- 
der the same assumptions as in the Eq. (|15fl is given by 
the expression 



F = 1- 



32a, 

T7 



^| C o| 2 | Cl | 2 -^( C 2 + c?) (cf + cf) 



(21) 

From this expression we can conclude that when the ini- 
tial state vector is proportional either to |0) or |1) the 
contribution of the control errors made in the (y, r\) \e 1= o 
plane can be neglected at the accepted approximation de- 
gree. In our previous work [2l| the fidelity was defined as 
/ = \f~F~ and the condition |co| 2 |ci| 2 = was assumed. 
In the limit (r^^) -1 — * 0, when the fluctuations average 
out, from the Eq. (|2H we obtain that 1 — / ~ a x . Thus 
our previous result |2j| concerning the cancellation of the 
squeezing control errors up to the fourth order on their 
magnitude is reproduced as it should be (remind that a x 
has the order of (5r x ) 2 ). 

We can simplify the expression l|21l) exploiting the 
parametrization (|18fl for the amplitudes cq and c\ and 
assuming £ — \ = ™, where n is an integer. Under these 
assumptions the gate fidelity is given by 



F = 1 



8cr a 



T x 



^-T y sin 2 2<p. 
1 v 



(22) 



It is evident that the gate fidelity oscillates depending 
on the choice of the initial qubit state. It is less for 
the superimposed initial states than for the basis ones. 
Namely, the fidelity has its maximum when the initial 
qubit state |^o) is proportional either to |0) or |1). In this 
case the errors made in the (y,r\) \e 1= o plane are elim- 
inated. As well the fidelity has its minimum when the 
initial state vector is equal to (|0) ± |1)) /V2. Moreover 
from Eqs. I|17|l and (|21|l we find a simple linear formu- 



lae connecting the purity of the final state and the gate 
fidelity: 



2F- 1. 



(23) 



This expression demonstrates the close connection be- 
tween these quantities defining the decoherence strength 
and the gate stability. 



V. CONCLUSIONS 

We considered optical and ion trap HQC proposed the 
Refs. Q and respectively. Regarding the particular 
implementation of Hadamard gate we have studied de- 
coherence induced by stochastic squeezing control errors. 
Ornstein-Uhlenbeck stochastic process was exploited to 
model random fluctuations of the squeezing control pa- 
rameter. We have analytically obtained the purity of the 
final qubit state and calculated the fidelity of the non- 
ideal Hadamard gate. It was shown that the stochastic 
squeezing control errors reduce the final state into a mix- 
ture of pure states and, thus, induce decoherence. We 
have shown that the final state purity oscillates depend- 
ing on the choice of the initial qubit state. The purity 
decreases when the initial qubit state is a superposition 
of the basis states. For example, it has its minimum 
when the initial state vector is equal to (|0) ± |1)) /V2- 
In contrary, the purity has its maximum when the ini- 
tial qubit state is proportional to one of the basis states 
|0) or |1). In this case the control errors made in the 
(y,ri) |e 1= o plane are eliminated. The same conclusions 
can be made for the gate fidelity. Simple linear formulae 
connecting the gate fidelity and the purity of the final 
state was derived. 
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